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We solve the two-channel Anderson impurity model using the Bethe-Ansatz. We determine the 
ground state and derive the thermodynamics, obtaining the impurity entropy and specific heat over 
the full range of temperature. We show that the low temperature physics is given by a line of fixed 
points decribing a two-channel non Fermi liquid behavior in the integral valence regime associated 
with moment formation as well as in the mixed valence regime where no moment forms. We discuss 
relevance for the theory of UBei3 . 
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In recent years a large number of alloys have been 
observed to deviate from the Fermi liquid behavior, 
with the low temperature thermodynamics and trans- 
port physics described by logarithmic or fractional power 
laws. Among these are heavy fermion materials based 
on Ce 3+ or U 4+ ions containing inner shell /-electrons 
that do not delocalizejlj]. We shall concentrate hence- 
forth on U based materials, in particular UBei3. Hund's 
rules and spin-orbit coupling in the presence of a cubic 
crystalline electric field lead to modeling of the U ion 
by a Kramers doublet in a 5/ 3 configuration to be 
represented by the fermionic creation operator ft and 
a quadrupolar (non-magnetic) doublet T3 in 5/ 2 config- 
uration represented by the operator b\. The doublets 
hybridize with conduction electrons in r 8 representation 
carrying both spin (er =t, I) and quadrupolar (a = ±1) 
quantum numbers. Strong Coulomb repulsion requires 
single occupancy of the localized levels, ft fa + b\b & = 1. 
The resulting hamiltonian is the two channel Anderson 
impurity model©, |, |,|,|, 
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The spectrum is linearized around the Fermi level and 
the Fermi velocity is set to one with the resulting den- 
sity of states being p = 1/ (2tt) . We shall study the 
model in the grand canonical ensemble with the chem- 
ical potential fi coupled to the total number of electrons, 
N = f ip^ipav dx+ftf a . The magnetic and quadrupolar 
impurity doublets are at energies e s and e q respectively. 
The last term gives the hybridization of the host with 
the impurity. The bar over a indicates that the index 
transforms according to the conjugate representation. 

The model has been intensely studied recently via a 
new Monte Carlo method Q and by conserving slave bo- 
son theory Q. However, many open questions remain. In 



this letter we shall present the solution of the model via 
a Bcthc-Ansatz construction, applicable in the past only 
to the single channel case f|. We shall give a complete 
determination of the energy spectrum and the thermo- 
dynamics, allowing us to follow the evolution of the im- 
purity from its high temperature behavior with all four 
impurity states being equally populated down to the low 
energy dynamics characterized by a line of fixed point 
hamiltonians H*(e, A) where e = e s — e q and A = irpV 2 
(we shall hold A fixed in what follows). 

We shall find that the line of fixed points is charac- 
terized by a zero-temperature entropy Sf mp = fcsln\/2 
and a specific heat C™ p ~ T In T typical of the 2-channcl 
Kondo fixed point. However the physics along the line 
varies with e. Consider n c = (ft fa), the amount of 
charge localized at the impurity. For e < /j — A, we find 
n c k, 1 signaling the magnetic integral valence regime. 
At intermediate temperatures a magnetic moment forms 
which undergoes frustrated screening as the temperature 
is lowered, leading to zero-temperature anomalous en- 
tropy and anomalous specific heat. For e > p+ A it is the 
quadrupolar integral regime and a quadrupolar moment 
forms. In the mixed valence regime, |e — p\ < A, similar 
low temperature behavior is observed though without the 
intermediate regime of moment formation. In more de- 
tail, each point on the line of fixed points is characterized 
by two energy scales, T;^(e). These scales describe the 
quenching of the entropy as the temperature is lowered: 
the first stage taking place at the high temperature scale 
Th, quenching the entropy from fc^ln4 to feeing, the 
second stage at 7], quenching it from ku In 2 to ks In 
In the integral valence regime, \e— p\ 3> A, the two scales 
are well separated and as long as the temperature falls 
between these values a moment is present, (magnetic or 
quadrupolar depending on the sign of e — /1) , manifested 
by a finite temperature plateau Si mp = ks In 2 in the en- 
tropy. It is quenched when the temperature is lowered 
below Ti, with T; — > Tk in this regime. For e = \x the 
scales are equal, Tj(/i) = Th{p), and the quenching oc- 
curs in a single stage. In a subsequent paper Q we shall 
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present the line of boundary conformal field theories cor- 
responding to our Bethe-Ansatz solution, and show that 
the two scales parameterize the approach to the fixed 
points. 

The Bethe-Ansatz wave functions consist of plane 
waves with momenta {%}, and amplitudes that are con- 
nected by a set of scattering matrices (S-matrices) de- 
rived from the Hamiltonian and obeying Yang-Baxter 
consistency conditions. The electron-impurity scattering 
matrix is, 
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where the index zero denotes the impurity. The op- 
erator Q acts on quadrupolar (or flavor) space as an 

"annihilation-creation" operator, [Q]° ^ = S a .jjS a . 
This S-matrix is unitary and group invariant. To find 
the electron-electron S-matrix we solve the two particle 
problem and find a matrix equation constraining the pos- 
sible form of the S-matrix, 

(82,0 — I)(Si,oSi,2 — I) = P? ,2-^1,2(^1,0 — I)(S2,0 — Sl,2) 

where P 9 is the flavor exchange operator and V s the spin 
exchange operator, with [P 9 ]^.^ = 8@ 5p and similarly 
for P s . Imposing in addition the Yang-Baxter conditions, 
we obtain an overdetermined set of equations admitting 
nevertheless a unique solution: 
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(fci - k 2 ) - i2A Pf _ 2 {h -k 2 ) + i2A P« 2 
(fci - k 2 ) - i2A (fci - k 2 ) + i2A 



The Yang-Baxter conditions and the single occupancy 
constraint allow the consistent generalization of the so- 
lution to an arbitrary number of electrons, N. Impos- 
ing periodic boundary conditions, one is led via standard 
methods to Bethe-Ansatz equations (BAE) that allow the 
determination of the momenta and hence the spectrum: 
E = T,f=ikj- The BAE are 
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where e„ (z) = (z — in A) / (z + in A). M s is the number 
of down spins and the spin-rapidities Af, describe their 
dynamics. The same role is played by M q and the fla- 
vor rapidities A 9 . The momenta {kj} (charge-rapidities) 
analogously describe the charge dynamics. 

Analyzing the equations in the thermodynamic limit 
we find that the solutions fall into four classes: (i) real 



charge-rapidities; (ii) charge-spin strings: k± = A ' s ±iA; 
(iii) spin-strings: A^ (ms) = A a ( m ") + iA (1 + m s - 2n), 
n = 1,... , m s (with arbitrary length m s ); (iv) flavor- 
strings: A^ (m,) = A a( - m «1 + iA(l + m g -2n), n = 
1,... ,m q (with arbitrary length m q ). A generic full- 
solution consists of a combination of all four kinds. 

Identifying the ground state and the excitations, we 
sum over the latter to derive the free energy. It 
is expressed in terms of an infinite set of functions 
of the variable A = j^(k — /1). These functions, 
{%(A),r? u (A),r?"(A),r?™(A)} n= i i ... ;00 , determine the dis- 
tribution probabilities of the various excitations at tem- 
perature T (entering the equations as r = the 
charge bound states, charge unbound states, spin excita- 
tions and flavor excitations, respectively. These distribu- 
tion functions satisfy an infinite set of coupled non-linear 
integral equations: the thermodynamic Bethc Ansatz 
(TBA) equations, 
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At 

lnry s „ = G * [S n ,iU + <Wfsn-i + isn+i] 

lnffan = G * [—Sn t ifn - 5 n ^ + <5„,ifg„_i + f gn +i] 



where we have introduced the notations 77 = 1/r], i a = 
In (I + rja), f a = ln(l + fj a ) and <5„,i = 1 - 5 n ,i- The 
convolutions (denoted by *) are defined with the ker- 
nel G(A) = 2 - c ^ sh A . The magnetic field h s and the 
quadrupolar field h q determine the asymptotic condi- 
tions: lim„^ 00 (.K„ + i * f s „ - K n * f sn+ i) = -2h s /T and 



lim„^oo(i^ n+ i * t qn - K n * iqn+l) 



-2h q /T, where 



K n (\) — ( 2A )2 4 "(„ 7r )2 • 

In terms of the distribution functions, the free energy 
F = -Fbuik + ^imp is given by: 

F bulk = —TL I pfudX- 2TL I pi- b d\ 



imp 



e q - T [G * (K x * f„ + K 2 * f 6 + f ?1 )] ; 



where J = Having solved for the distributions for 

given values of temperature and fields, one can compute 
the free energy of the impurity for any value of e. The 
bulk part of the free energy is divergent and requires the 
introduction of a regularization scheme. On the other 
hand, the impurity contribution to the free energy is reg- 
ular and independent of the cut-off procedure. 

We turn now to study the TBA equations, considering 
first the zero temperature limit for which it is convenient 
to introduce £ = Tin n. We then have, lim r _> + ri = £ + 
and lim r ^ + rf = £~ , where 2£ ± = £ ± in terms 
of which we are able to solve explicitly the TBA equa- 
tions. We find: £ fc = A, £ q2 = £~ 2 = G * ^ , = 
= -| — £ q2 with all the other £'s vanishing (the re- 
sult of the convolution giving t; q2 can be written down 
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in terms of dilogarithms). Thus in the absence of ex- 
ternal fields, the ground state is built out of a sea of 
charge-spin strings filled up to A = (i.e. k = fi) 
and a completely filled sea of 2-flavor-strings. The zero- 
temperature impurity level occupancy (and therefore also 
the charge susceptibility) can be deduced in a closed 



form: 



= r 
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The occupancy is 



r+oo i(X-7r/J)j q2 , x 
[(A— 7r/ J)'- 1 +ir' 2 }' 2 dA 

integral: ~ 1,0 for \e — ju| 3> A, and non-integral 
elsewhere. 

We turn next to study the low temperature physics, 
r <C 1. In the integral valence regimes, \e — fi\ ^s> A, 
\J\ 3> 1 and the main contribution to the free en- 
ergy comes from |A|«7r/|J|3>l. It is convenient to 
rewrite some equations in the TBA. Using the exact re- 
lation: 111% — Infjgi = A/r, we can eliminate fj from 
the equations for fj u and rj q 2 and write them as: hifj u = 
-Ef/r+G*[f sl + i ql ] and In ^ 2 = -E d _/r+G*[l ql + f q3 ] 

with E%lr = G*ln(l + e ±T( A -^ ln ^i)^ . The terms E% 

become driving (inhomogeneous) terms at low tempera- 
tures as rln^gi tends to zero (with r 2 corrections). We 
can further approximate them as follows: 
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Changing variables, £ = A— j, the driving terms become, 
in the integral valence regime: E±/t = e ± (- inT / T ± 
where T± = ±0-e ±7T / J . The relation of T± to T Ul de- 
pends on the sign of e — /i and is discussed below. 

We analyze separately the magnetic and quadrupolar 
moment regimes. Consider first the magnetic case where 
e <C M — A, and therefore A <C — 1. In the limit of small 
t, the driving term E_ diverges faster than E+ that is 
compensated by a decaying exponential in the numera- 
tor; thus rj q 2 tends to zero exponentially fast cutting away 
the equations for the higher q-flavor 77's. After the identi- 
fications fj q i — > rjf, fj u — ► r)2 and rj sn — ► ?/* +2 , we recover 
the TBA equations of the 2-channel Kondo problem ||. 
The Kondo temperature is Tk =Ti = T + . (T_ is large, 
outside the low temperature approximation range). The 
identification of the resulting TBA equations in this limit 
with those of the 2-channel Kondo model indicates that 
at low temperatures, when E'L/t becomes very large, the 
system has a localized magnetic moment. When temper- 
ature goes below Tk, overscreening will take place and 
in particular an entropy Simp = ks ln-s/2 will arise. This 
will also be found numerically (see below). 

On the other hand, when e ^S> /i + A, the limit of 
vanishing r drives fj u to zero exponentially fast and cuts 
away the spin 77's. No remapping is required, and we re- 
cover again the TBA equations of the 2-channel Kondo 
model but this time for the quadrupolar degrees of free- 
dom. The relative minus sign in the driving term reverses 
the 'direction of lowering temperatures' with respect to 
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FIG. 1: Impurity contribution to the entropy as a function of 
temperature for different values of e — /1. As the temperature 
goes to zero all curves approach the universal value ks In \[2. 
Positive and negative values of e — /1 fall on top of each other. 



the previous case and has no further consequences. The 
Kondo temperature is Tk = T; = T_, below which the 
two spin channels will overscreen a localized quadrupolar 
moment. 

Now we study the mixed valence regime, |e — /i| < A, 
where n c w 1/2. As \J\ 1, we need the solutions 
around A ~ 0. Then both driving terms are approxi- 
mately E±/t = e ±C-inT/T ± anc j kp^jj diverge simultane- 
ously in the low temperature limit. We repeat the stan- 
dard procedure of shifting variables with ±lnT/TR-[||. 
Depending on the choice of sign, one of the driving terms 
diverges as 1/r 2 and the temperature 'disappears' from 
the other one. We recover one of the two cases out- 
lined above but this time the equations are valid only at 
temperatures below T/ and there is no localized moment 
regime. However, a residual entropy S^mp = ks 

In ^2 

arises also in this regime. 

Finally, we turn to the finite temperature thermody- 
namics. We obtain it by solving numerically the TBA 
equations. In Fig. [l] we show the behavior of the im- 
purity entropy as a function of temperature for different 
values of e. At high temperatures the entropy is ks In 4 
in agreement with the size of the impurity Hilbert space. 
For \e — /i| 3> A the impurity entropy is quenched in 
two stages. The degrees of freedom corresponding to the 
higher energy doublet are frozen first. The entropy be- 
comes fcsln2 and the system is in a localized magnetic 
or quadrupolar moment regime depending on the sign of 
e — [i. As the temperature is further decreased, the re- 
maining degrees of freedom undergo frustrated screening 
leading to entropy ks In V%- On the other hand, for val- 
ues of |e — A, the quenching process takes place in a 
single stage. As e — /1 is varied the behavior interpolates 
continuously between the magnetic and the quadrupo- 
lar scenarios. The zero temperature entropy is found 
to be independent of e in accordance with the analytic 
study of the TBA equations. Note also the presence of 
a crossing point, a temperature T cross w 0.1A where all 
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FIG. 2: Main plot: impurity contribution to the specific heat 
as a function of temperature for different values of e — fi. As 
this parameter approaches zero the Kondo contribution (left) 
and the Schottky anomaly (right) collapse in a single peak. 
Inset: experimental data for the 5/-derived specific heat of 
UBei3. 



entropies take the value Simp = kg In 2 independently of 
e (cf. Ref. P). 

In Fig. Jwe show the impurity contribution to the 
specific heat. The two stage quenching process gives rise 
to two distinct peaks. The lower temperature peak is 
the Kondo contribution centered around Ti whereas the 
higher temperature peak often referred to as the Schottky 
anomaly is centered around a temperature Th- Approxi- 
mate expressions for Thj can be read off from the curves: 
T h ,i{e) w ^jln(l + 2ae ± ^ e ~^) with 1 < a < 4 (this 
expression goes over to T± defined before in the appro- 
priate limits). For large \e — [x\ the two peaks are clearly 
separated and the area under the Kondo peak is ks In y/2 
while that under the Schottky peak is ks In 2. 

As mentioned earlier, the model was proposed as a de- 
scription for the U ion physics of UBei 3 . It is expected to 
describe the lattice above some coherence temperature. 
We provide in the inset of Fig. || the experimental data 
for the 5/-derived specific heat of the compound. It is 
obtained by subtracting from its total specific heat, the 
specific heat of the isostructural compound ThBei3 con- 
taining no 5/ electrons fl2"|. This way one is subtracting 
the phonon contribution as well as the electronic con- 
tribution from electrons in s, p and d shells (the pro- 
cedure is quite involved and we refer the reader to the 
articles cited just above for full details). The sharp fea- 
ture at 0.8K signals the superconducting transition of 
UBei3 and falls outside the range where this compound 
might be described by an impurity model. 

Concentrating on the temperature range containing 
the Kondo and Schottky peaks we conclude that no val- 
ues of e and A yield a good fit. Further, the entropy 
obtained by integrating the weight under the experimen- 
tal curve falls between fee In 4 and fcs In 6. This suggests 
that a full description of the impurity may involve an- 



other high energy multiplet (possibly a triplet cf. Ref. ||), 
almost degenerate with the T3 to yield a single peak 
for the Schottky anomaly. The nature of the multiplet 
could be deduced from further specific heat measure- 
ments. For an n-plet degenerate with the T3, one has 
an SU (2) x SU{n + 2) Anderson model and the area un- 
der C™p/T is then given by k B m[^±A sec while if 
the n-plet is slightly split off the doublet, the area is given 
by ks M 1 ^ 1 ] 0- In order to be certain of having elim- 
inated lattice effects it would be better to carry out the 
measurements on Ui_ a; Th ;r Bei3. For x > 0.1 the com- 
pound has no longer a superconducting transition and the 
lattice coherence effects are largely suppressed. Further, 
there are several experimental indications that support 
the idea of an impurity model description of the thori- 
ated compound for a wide range of temperatures (T^| . 

In subsequent work we shall present the solution of the 
general SU (JV) x SU (M) model and study the effects of 
magnetic and quadrupolar fields. These will help identify 
what particular impurity model is best fit for describing 
Ui-xThxBeia @- 
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